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Abstract 

C Consider the flow of a thin layer of non-Newtonian fluid over 
soUd surface. I model the case of a viscosity that depends 



|nonlinearly on the shear-rate; power law fluids are an impor- 
tant example, but the analysis here is for general nonlinear 
j^diependence. The modelling allows for large changes in film 
^^thickness provided the changes occur over a large enough 
^^lateral length scale. Modifying the surface boundary con- 
C^dition for tangential stress forms an accessible base for the 
a nalvsis where flow with constant shear is a neutral critical 
^^mode, in addition to a mode representing conservation of 
^^fluid. Perturbatively removing the modification then con- 
•structs a model for the coupled dynamics of the fluid depth 
'"^and the lateral momentum. For example, the results model 
C^the dynamics of gravity currents of non-Newtonian fluids 
^eyen when the flow is not very slow. 

,__|Keywords: thin fluid flow; non-Newtonian fluid; inertia; 
^power law rheology 
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^1 Introduction 

O • . 

■^Consider the two dimensional flow of a thin layer of fluid 
'^over a flat substrate. The fluid of thickness ri(x, t) spreads 
^with mean lateral velocity ii(x,t). Suppose the fluid has 
^tlie non-Newtonian, power law, stress-strain relation that 
*the stress oc (strain-rate)'' for some fixed exponent s: the 
exponent s — 1 for a Newtonian fluid; s < 1 is shear thinning; 
[and s > 1 is shear thickening. Such a power law is sometimes 
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hicalled Ostwald's or Norton's constitutive relation 5 . Then 



the systematic analysis developed in this article supports the 
nondimensional model 
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where Re is the nondimensional Reynolds number, Cs is the 
coefficient of proportionality in the nonlinear stress-strain re- 
lation, and where gi and gi are the nondimensional compo- 
nents of gravity along and normal to the flat substrate. Fluid 
is conserved through ([T]) . The momentum equation ([2]) incor- 
porates effects of inertia, self-advection, bed drag and grav- 
itational forcing; the dependence of the coefficients upon s 
models the subtle effects of the power law rheology. 

This model not only applies to the flow of simple liquids, 
it applies to: gravity currents of suspensions with medium to 
high volume fractions as these are non-Newtonian [17j ; power 
law rheologies are used to model ice flow O [THl e.g.] and at 
even a few metres per year the Reynolds number is signifi- 
cant for a thick glacier; and a modified model would apply 
to turbulent flow as the Smagorinsky large eddy closure of 
turbulence corresponds to the shear thickening case of expo- 
nent s = 2 Pl Eqn. (6), e.g.]. This article puts models such 
as ([I])-© within the sound support of modern dynamical 
systems theory. Section [H to empower us to systematically 
control error, assess domains of validity, and to systemati- 
cally account for further physical effects. 

The analysis here encompasses not only power law fluids 
but a general nonlinear dependence of the stress upon the 
strain-rate as codified in Section [H In contrast, almost all 
previous thin fluid film modelling use only a power law de- 
pendence. Some industrial plastics have a complicated non- 
monotonic dependence ^ that cannot be represented by a 
simple power law. Similarly, dense suspensions often have 
non-monotonic dependence [17j . The resultant model de- 
rived in Section |4] also applies to such complicated industrial 
plastics and dense suspensions. 

The lubrication approximation of very slow flow, low 
Reynolds number, underpins previous theoretical models for 
non-Newtonian thin fluid films: Perazzo & Gratton [T^ and 
Betelu & Fontelos [T] examined flow with surface tension; 
this followed experiments compared with travelling waves 
and similarity solutions by Gratton, Minotti & Mahajan [5]. 
Gratton et al. comment "the differences between Newtonian 
and non-Newtonian currents are significant and can clearly 
be observed in experiments" . But the lubrication approxi- 
mation, that creates models expressed only in terms of the 
fluid thickness T)(x,t), does not model inertia and so can- 
not resolve any wave-like dynamics. To model faster fiows, 
potentially with wave effects, we must resolve the dynam- 
ics of both the fiuid thickness and a measure of horizontal 
momentum |12[ [TC], we used rj and ii in ([I])-®. For exam- 
ple, Harris et al. [6^ modelled particle driven gravity currents 
using shallow water equations that resolve the dynamics of 
both the fluid thickness and the mean lateral velocity. How- 
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ever, such modelling of essentially dissipative flows, albeit 
dissipative via turbulence, by the laminar inviscid founda- 
tion of shallow water equations appears a contradiction that 
demands resolution. This article shows how such models of 
non- Newtonian fluid flow may be put on a sound mathemat- 
ical basis to empower accurate physical forecasts. 

2 Differential equations to model 
non-Newtonian flow 

Let the incompressible fluid have thickness r\[x,t), constant 
density p, a nonlinear rheology, and let the fluid flow with 
some varying velocity fleld u — (u,v) — (ui ,U2) and pres- 
sure field p. 



Nonlinear constitutive relation 
tensor [gl fTTpl 



Define the strain-rate 



9ui 
9xi 



9uj 
9xi 



(3) 



where xi = x and X2 — y are distances along and normal 
to the solid substrate, respectively. Then the stress tensor 
for the fluid is ffij = — p6ij -|- 2p"V£ij : for a Newtonian fluid 
the kinematic viscosity y is constant; but when the kine- 
matic viscosity varies with strain-rate then we model shear 
thickening or shear thinning non-Newtonian fluids. 

The important class of non-Newtonian fluids that we ad- 
dress has viscosity which depends only upon the magnitude £ 
of the second invariant of the strain-rate tensor [1] : 
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For example. Bird et al. [21 see [T]] report that a solution of 
0.5% Hydroxyethylcellulose is shear thinning: at 20 °C the 
solution has viscosity [i = me^^^ for exponent s — 1/1.96 
and coefficient m = 0.84 N s*^ / m^ . 

Partial differential equations Make equations nondi- 
mensional with respect to some velocity scale, a typical fluid 
thickness, and the fluid density. The nondimensional pdes 
for the incompressible, two dimensional, fluid flow are firstly 
the continuity equation 
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and secondly the momentum equation 
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where Re is the appropriate Reynolds number, T is the 
nondimensional deviatoric stress tensor, and g — (91,92) 
is the nondimensional forcing of gravity. For a fluid with 
a nonlinear stress-strain relation, the nondimensional devia- 
toric stress tensor 
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^Some, such as Betelu &; Fontelos [T], use double this tensor. 



Boundary conditions Solve these pdes with nondimen- 
sional boundary conditions: 



on the bed of no-slipH 



u 



on y = ; 



the kinematic condition on the free-surface of 
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• the stress normal to the free surface comes from con- 
stant environmental pressure and surface tension, that 
is, 
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where We is a nondimensional Weber number charac- 
terising the importance of surface tension; 

• and there must be no tangential stress at the free sur- 
face, 

(1 -ti^)ti2 -Ft1x(t22 -Til) =0 on y=ri. (11) 

This boundary condition of zero tangential stress im- 
plicitly is effectively one of zero shear at the surface; 
this is not appropriate for material with a flnite yield 
stress. Here we assume the fluid yields for arbitrarily 
small stress. 



3 Centre manifold theory supports 
the modelling 

This section shows how to place models such as (HJ-© on 
a sound theoretical base. Artiflcially modify the zero tan- 
gential stress free surface condition (fTT|) to have an artiflcial 
forcing proportional to the local velocity, a forcing which we 
later remove by evaluating at parameter y — ^ : 

(1 - It) [d -T1^)Tl2+T1x(T22-Tii)] 

= (l-y)— on y=ri. (12) 

Evaluated at y = 1 this artiflcial right-hand side becomes 
zero so the boundary condition (|12p reduces to the physical 
boundary condition (jlip . However, when the parameter y = 
and the lateral gravity and lateral derivatives negligible, 
91 = 9x = , a neutral mode of the dynamics is the lateral 
shear u — v^Ey where I define E to be proportional to the 
mean lateral strain-rate: 
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^If modcUing turbulent flows by a large eddy closure, we may justi- 
fiably replace this no-slip bed condition by a mixed boundary condition 
on the lateral velocity: ti oc . 
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This neutral lateral shear mode arises because in pure shear 
flow Ti2 = yuy and hence the artificial free surface condi- 
tion reduces to "vUy — yu/r\ on y — r\ . Conservation 
of fluid provides a second neutral mode in the dynamics. 
That is, when Y = gi = 9x = then a two parameter 
family of equilibria exists corresponding to some uniform 
lateral shear flow, u = Ey , on a fluid of any constant thick- 
ness r\. For large enough lateral length scales, these equilib- 
ria occur independently at each location x [TTl [131 e-g-] and 
hence the space of equilibria are in effect parametrised by 
e(x] and ri(x). Provided we can treat lateral derivatives 9x as 
a modifying influence, that is provided solutions vary slowly 
enough in x, centre manifold theory [4j[7l[T4j e.g.] assures us 
three things: this space of equilibria is perturbed to a slow 
manifold, on which the evolution is slow, that exists for a 
finite range of y and gi, and which may be parametrised 
by the mean lateral shear E(x,t) and the local thickness 
of the fluid ri(x,t); the slow manifold is attractive for all 
nearby initial conditions; and that a formal power series in 
the parameters y, gi and derivatives 9^ approximates the 
slow manifold. That is, the theory supports the existence, 
relevance and construction of slow manifold models such as 



4 Low order models of the dynamics 

Computer algebra readily constructs slow manifold models 
as asymptotic solutions of the governing differential equa- 
tions and boundary conditions. 



ity fields computed simultaneously with (fT5 |) -(fn | the com- 
puter algebra [15, §3] also derives the mean lateral velocity 
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Reverting this series to express e in terms of ii, and sub- 
stituting into the model p ^ -lfTl ]) leads to a model for the 
coupled evolution of ri(x, t) and u(x, t). Evaluating at the 
physically relevant y = 1 , to remove the artifice in the sur- 
face boundary condition ((T^), then gives the model II])"® 
discussed in the Introduction of this article. 

Computer algebra experiments [151 §1] suggest that the 
convergence of the asymptotic series in y is markedly im- 
proved by the factor (1 — -^y) on the left-hand side of the 



tangential stress boundary condition ([T2|) . This factor is 
equivalent to an Euler transformation of the asymptotic se- 
ries. As shown in other similar applications [T^ e.g.], eval- 
uation at Y = 1 is valid provided the lateral derivatives are 
small enough. 

Computer algebra |15l §3] may construct terms in the for- 
mal power series solutions to higher order in the parameters 
Y, gi and 9x to generate many valid approximations of vary- 
ing orders of accuracy. For example, to resolve any effects 
of surface tension we need to compute terms in 9^ that are 
neglected in ^ and With the support of centre man- 

ifold theory, researchers may choose an approximate model 
that suits the parameter regime of their application. 



4.1 Power law fluids 

For simplicity, suppose the rheology is a nondimensional 
power law for the kinematic viscosity, y = c^t^^^ ■ 

Computer algebra [15, §3] derives that for such a power 
law ffuid, the evolution of the ffuid thickness r| and the stress 
parameter E is 
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The nonlinear rheology primarily appears as a nonlinear drag 
on the bed. However, changes to the vertical profiles of ve- 
locity and pressure due to different power laws affect the 
coefficients of this model through their dependence upon ex- 
ponent s. 

In modelling the ffow of thin ffuid layers, researchers gen- 
erally prefer to use the mean lateral velocity or the lateral 
ffuid flux instead of the shear parameter E. Using the veloc- 



4.2 More general non-Newtonian fluids 

We now return to the more general case where the viscosity "v 
of the fluid depends quite generally upon the magnitude of 
the shear-rate £, instead of being a simple power law. In 
this more general case the expressions for the modelling are 
much more complicated. The reason is the general nonlinear 
dependence of viscosity on strain-rate: for conciseness deffne 
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where primes on 'v denote the derivatives d/ dE of the vis- 
cosity -v(e) and evaluated at e = \/2u/r\ . 

The procedure is as for the power law case: computer al- 
gebra [IHl §3] constructs the slow manifold and evolution 
thereon to some order of error; then revert the asymptotic 
series to find stress parameter E as a function of mean ve- 
locity u; and substitute to express the model in terms of 
ri and u. Conservation of fluid again derives ([T]) (to any 
order of error). The dynamics of momentum then leads to 
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Evaluate this equation at y = 1 to recover a physically rele- 
vant model of the dynamics of lateral momentum. 

The power law model ^ is just one specific subclass of 
the general model (fTn|) : obtain ^ by the specific choice of 
a power law viscosity -^(e) — Csi^^^ ■ 

5 Conclusion 

Following similar modelling for Newtonian thin films [12] . 
this approach places the modelling of an important class of 
non-Newtonian fluids upon the sound basis of centre mani- 
fold theory [U [3 [TH e.g.]. This modern dynamical system 
foundation empowers us to systematically derive the novel 
and accurate models ((Sj), (fT4)) and (fT6)) for the lateral mo- 
mentum of fluids with nonlinear rheology. 

These models of thin fluid flow can be directly applied 
to flows as diverse as those of industrial plastics ^ e.g.], 
ice [H1[TS1 e.g.], and medium to dense suspensions [TTl e.g.]. 
When you desire more accuracy than that presented here, 
computer algebra readily computes higher order approxima- 
tions [TS] §3]. Modifying the no-slip boundary condition on 
the bed, ([5]), will empower the modelling of turbulent lay- 
ers of flow over a substrate via large eddy closures. There 
are enormous applications for this approach to modelling the 
dynamics of laterally extensive layers of fluids. 
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